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Summary: - We first obtain by analogy with the continuous (differential) case the gen- 
eral solution of the discrete Riccati equation. Moreover, we establish the full equivalence 
between our discrete Riccati equation and a corresponding homogeneous second order 
discrete linear equation. Our results can be considered the discrete analog of Miclnik's 
construction in supersymmetric quantum mechanics [J. Math. Phys. 25, 3387 (1984)]. 
We present an application to the three-site master equation obtaining explicitly the 
general solutions for the simple cases of free random walk and the biased random walk. 

PACS 11.30.Pb - Supersymmetry 

We consider the continuous Riccati equation (CRE) y — a(x)y 2 +b(x)y + c(x) with 
the known particular solution ?/o and let y\ = u + yo be the second solution. By 
substituting yi in CRE one gets the Bernoulli equation u = au 2 + (2ay + b)u. 
Furthermore, using v = l/u, the simple first-order linear differential equation v + 
(2ayo + b)v + a = is obtained, which can be solved by employing the integration 

factor / = e r 2ayo+b , leading to the solutions = e ~ f 2ayo+b (- f aef 2aVo+h +C) , 
where C is an integration constant. Thus, the general solution of CRE is 

e J 2ayo+b j 

vi = va + pr- — - = yo + — — r x, ■ (1) 

For the usual nonrelativistic one-dimensional quantum mechanics 6 = and a = — 1, 

i.e., the integration factor is f qm = e~ 2 / Vo . The particular solution yo is known as 

Witten's superpotential jjj in the context of supersymmetric quantum mechanics Q , 

while the general solution y% = y$-\ fe may be called Mielnik's superpotential 

c+J f qm 

||. Moreover, f^m is the ground state wavefunction, — 2y (= —2^3 Inf^m) is the 

/ 1/2 

Darboux transform part of the Schrodinger potential and — HH? is the ground 

c+J f qm 

state wavefunction corresponding to Mielnik's superpotential for which —2y 1 can be 
thought of as the general Darboux transform part in the potential. All these rela- 
tionships make up the core of the "entanglement" between Riccati and Schrodinger 
equations, which has been recently emphasized by Haley [Q. 

Consider now the discrete Riccati equation (DRE) y n +i = a n yn+iyn + b n yn + Cn, 
with known solution As in the case of CRE, the general solution is chosen of the 
form y n = y n + u n . The discrete equation for the function u n (discrete Bernoulli 



equation) will be 



a n b n + a n y n +i 
U n +1 = TrU n+1 U n H — u n . (2) 



We now take v n = 1/ u n to get 



( 1 - a nV n \ a n 
Vn + 1 = I 1 Vn ~ I ; = SnV n + K , (3) 



whose solution reads 



k 



-l 



\fc=0 /V fe=0 \j=0 ) ) 

where D is a positive constant. Thus, the general DRE solution can be written 
down as follows 

T-fn-l ^ b k +a k y k + 1 . 

n^t^^) M 

_ „_i / fe fej+Ojj/° \ q fc 
^ Z^fe=olllj=0 l-a i2/ » jb k +a k yl +1 

One can remark the formal analogy between the CRE and DRE general solutions. 
Indeed, f p = IlzLo b '^ lVl ^ 1 can be thought of as a discrete "integration" factor 
(DIF), and thus Eq. (5a) can be written in the form 



, J _ , ,0 _| fn 

ri _ y^™ -1 f k a k 

^k=0 b k +a k y» +1 



Vn = Vn+ ~ 7—— , (56) 



from which the analogy (and difference) with respect to the continuous case is 
transparent. 

We now want to establish a formal equivalence between the DRE and the homo- 
geneous second-order discrete equation (HSODE), more exactly, between y n +i = 
anj/n+ij/n + b n Vn + c n and c' n x n+2 + b' n x n+ i + a n x n = 0. 

In the DRE let y n = -^-^ + d n . After some easy algebra and with the crucial 
choice d n+ \ = —b n /a n one gets 

(c„ + —)x n+2 - (1 + a n ^— -)x n+ i + a n x n = . (6) 

^n— 1 

Thus, the identifications are b n = —1 — and c n = c n + ^ £L . Vice versa, let us 

start with the HSODE and introduce y n = x n /x n+ i + (b n + l)/a n . Again, after quite 
simple algebra, one gets the DRE with the following identifications b n = —a n (l + 
b n+1 )/a n+ i and c n = c n + (1 + b n+1 )/a n+ \. These are the inverse identifications of 
the first ones. 
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One may foresee many physical applications of the above mathematical proce- 
dure. Here, we present only one application, namely to write down the general 
solution of the Hermitian three-site master equation written as in 

@ and @ 

~{gn+ir n +2) 1/2 x n+ 2 + (g n +i + r n +i)x n+1 - (g n r n +i) 1/2 Xn = , (7) 

where g n is the transition rate for the jump n — > n + 1 and r n is the one for the 
backward jump n — ► n — 1, while be thought of as the square root of the 

site probability. The known stationary solution is 



•>(1I^)' J - 



r 

. j=o ■ 



where A is a scaling constant, that through the normalization condition can be fixed 
to be A = f 1 + J2n =1 n™=i H^j ( see !)• The transformation 

_ X n _ 1 + gn+l + Tjg+l , , 

in+i {gnr n +i) ' 

applied to Eq. (7) leads to the following master Riccati equation 

/ \ 1/2 

( sl/2 / ffn 7 'n+l \ , . s ( sl/2 1 + .9n+2 + 

2/n+i = - .9nr, i+ i 1 yn+iy n ~\ (l+g n+ 2+r n+2 )y n -{g n+ ir n+ 2) 1 — --j— 

\gn+\r n +2 J {g n +ir n +2Y ,z 

(10) 

with the particular solution 

y0 = 1 ± g " +1 (11) 

an i M /2 ' \ xx ) 



However, one can write the general solution corresponding to Eqs. (5a, b), where 

xf ee -(g k r k+l )V* and 6 fe = bf 



y° is substituted by Eq. (11) and a k — a k ! = -{gtrk+i) 1 ^ 2 and b k = b^f ~ 



1/2 



(1 + gk+2 + ^fc+2), where the superscript M stands for master. By 



\ Sk+irk+2 

introducing the master DIF 

1 ^-^ 1 — a, vr 

1=0 1 yi 

one can write a compact form of the general master Riccati solution 

, M (9on) 1/2 , x° s 2 

1 _ , fn-l _ 1 + ffn+1 9nTn+l y xl + 1 > 

(13) 
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where D is a constant. Moreover, using simple discrete algebra, one can obtain the 
following new general solution of the master equation 



n-l / fegri^(_g£_)2 x-1 

4=4TT ! + ^ o , (14) 

where x° is the known product solution Eq. (8) having p factors. 

Let us see now the particular case g n = c\ , r n = c 2 , where c\ and c 2 are two real 
numbers G (0, 1), such that c\ < c 2 , i.e., q = ^ < 1. Then, the known stationary 

master solution is x" = Aq 1 ^ , the particular Riccati solution is = — (~^T72 , 
whereas the general Riccati solution reads 

■ ■i- 1 + ^ | g~" (15) 

(cic^D-g-V'E^oV* ' 

At this point, we treat separately the two subcases as follows. 

(i) For ci = c 2 = \ (free random walk), the particular solution is y° = —3 and 
the general solution is 

(ii) For ci ^ c 2 (biased random walk), say, c\ = |(1 — e) and c 2 = |(1 + e) so 
that q — j the particular solution is = — ^ 1 _ 3 ( ^^i/2 an d the general one is the 
following 



l 3-6 g-» 3-6 , (£§)" 



(1-62)V2 ^.^-l/a^p (1-6^)1/2 ^ + ( i ^)Vi [1 _(l±l)n] ■ 

£ (17) 

For e^O this solution goes into the above one for the free random walk. 

Passing to the master equation one can get in the two subcases the following 
solutions: 

(i) x" = A and x x n = A(l — jy^i)- For the positivity of the solution one should 
take D > n — 1. One can see that for D w n, at large n, x* is nought. However, 
for D = an, where a > 1, again at large n, x\ = A(l — — ) is substantially different 
from x°. Finally, for a 3> 1 one gets the particular solution. 



i-( 
d 



(ii) x° = A(^) "i 1 and after reabsorbing some constants one gets x\ — x° ^ ^ 

Inspection of the latter formula shows that for positivity one should take D > 



, where p = which is close to D > n + 1 resulting from the first 



- 1 , where p - 
order approximation. 

In summary, we derived the explicit form of the general solution of the three-site 
master equation based on the corresponding discrete Riccati general solution. In the 
continuous case it corresponds to Mielnik's procedure in supersymmetric quantum 
mechanics. Since the steady state solution is always a particular solution, what 



we have obtained here are either transient solutions H or solutions determined by 
some noncquilibrium boundary conditions fixing the value of the constant D. 
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